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Abstract 

Bouncing droplets on a vibrating fluid bath can exhibit wave-particle behavior, such 
as being propelled by interacting with its own wave field. These droplets seem to walk 
across the bath, and thus are dubbed walkers. Experiments have shown that walkers 
can exhibit exotic dynamical behavior indicative of chaos. While the integro-differential 
models developed for these systems agree well with the experiments, they are difficult 
to analyze mathematically. In recent years, simpler discrete dynamical models have 
been derived and studied numerically. The numerical simulations of these models show 
evidence of exotic dynamics such as period doubling bifurcations, Neimark-Sacker 
(N-S) bifurcations, and even chaos. For example, in [ ], based on simulations Gilet 
conjectured the existence of a supercritical N-S bifurcation as the damping factor in 
his one- dimensional path model. We prove Gilet’s conjecture and more; in fact, both 
supercritical and subcritical (N-S) bifurcations are produced by separately varying 
the damping factor and wave-particle coupling for all eigenmode shapes. Then we 
compare our theoretical results with some previous and new numerical simulations, 
and find complete qualitative agreement. Furthermore, evidence of chaos is shown by 
numerically studying a global bifurcation. 
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1 Introduction 

Pilot wave theory is based upon the idea that particle trajectories drive the statistics seen 
in quantum mechanics. In the early days of quantum mechanics it was seen as a promising 
explanation of the statistics arising in experiments. Specifically the works of de Broglie [ 
and later Bohm [3], while very different, showed potential as alternate formulations to that 
of the Copenhagen interpretation until the 1960s. In recent years, fluid dynamic experiments 
with droplets bouncing on a vibrating bath were observed to exhibit quantum-like behavior. 
This effectively revived the theory in the form of hydrodynamic quantum analogs and cat¬ 
alyzed the current very active state of walking droplet research. A summary of quantum-like 
behavior in experiments of walking droplets can be found in [ 1, 5, 6]. 
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In 2005 Couder and collaborators observed droplets moving across a vibrating fluid bath 
j ]. The droplets accomplish this by exciting an eigenmode at each bounce, which generates a 
wave field, and in turn are propelled by interacting with this field [7, 8]. Later, Couder et al. 

] and Eddi et al. [10, 11] observed quantum-like behavior such as single particle diffraction, 
tunneling, and the Zeeman effect. Furthermore, there were extensive experiments conducted 
on the quantization of orbits by Fort et al. [12], Harris et al. [13, 1 1], and Perrard et al. 
[15, 16], 

The experimental work led to many mathematical models. Discrete path memory models 
(accompanied by experimental results) for quantized circular orbits and walking in free space 
were developed by Fort et al. [ ] and Eddi et al. [ ], respectively. A detailed hydrodynamic 

model derived from free space experiments was developed by Molacek et al. [18]. Oza et al. 
took this further and developed integro-differential equation models for walkers in free space 
[19] and a circular rotating frame [20, 21], Then more accurate models for the waves were 
developed by Milewski et al. [22] and coupled with the equations of motion of the droplet. 

While many of these models agree very well with experiments, the equations are quite 
difficult to study analytically. The complexity of the equations naturally created interest in 
developing realistic simplified mathematical models exhibiting important dynamical features 
of the original, while being easier to analyze. The first models to exploit the discrete nature 
of the bouncing/walking were developed by Fort et al. [ ] and later Eddi et al. [ ]. Several 

other such reduced models were developed and two of them, devised by Shirokoff [23] and 
Gilet [1], showing considerable promise are actually planar discrete dynamical system models. 
Shirokoff [23] developed a model in which he derived a map for the motion of a particle in 
a square cavity. In this model, using numerics, he discovered cascading period doubling 
bifurcations indicative of chaos. Gilet [1] included the amplitude of subsequent modes in 
his model (1) for the straight line motion of a particle. He observed what appeared to be a 
Neimark- Sacker (N-S) bifurcation in numerous simulations, and so conjectured its existence 
and type. 


W n+ i = H[w n + T^)], 

X n+ 1 = X n - Cw n ty'(x n ); 

The model diagram is illustrated in Fig. 1. In this model w n G M is the amplitude of the 
wave just before impact n G N, x n G M is the position of the walker at impact n, C G [0,1] 
represents wave-particle coupling (related to the size of the droplet), G [0,1] is the damping 
factor (related to the path memory), and $ G 1 is a single eigenmode of the Faraday wave 
field. From the definition given in [1], one can calculate the damping factor of mode k with 
the associated memory Mk as )ik = exp(— 1 /Mk). However, since this model has only one 
eigenmode, k — 1. One may also interpret the damping factor as the ratio of the wave 
amplitude just before impact n + 1 to the amplitude just after impact n. Then, we think of 
fi = 0 as a system with no memory and hence no walking. Furthermore, fi = 1 corresponds 
to infinite memory and gives rise to Faraday waves; that is, there are waves present even 
when a drop is not. In addition, C = 0 corresponds to a drop that is nonexistent and C = 1 
represents a drop that is so large that the difference in position is exactly the gradient of the 
wave field at impact n. 
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Figure 1: Diagram of (1). 
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While it is not possible to construct an exact experiment for this model because con¬ 
straining the droplets to an annular or rectangular domain would cause boundary effects 
from the waves, we may think of this model as being analogous to or an approximation 
of some experimental setups. Two such experiments have been conducted by Filoux et al. 
24, 25]; the first being walkers in an annular domain and the second an upcoming work on 
walkers in a long-narrow rectangular domain. 

Although the focus of this investigation is rigorous mathematical analysis of certain 
aspects of the model dynamics, in order to connect the theory with the physics we shall 
discuss the physical implications of the dynamical results. Furthermore, the theory developed 
will give us insight into experiments such as [24, 25] and perhaps even 2-dimensional free 
space experiments due to the connections between discrete dynamical systems and continuous 
dynamical systems. 

The remainder of this paper is organized as follows: In section 2 we prove the basic 
dynamical properties of (1). We analyze the map via dynamical systems and bifurcation 
theory in section 3 and formulate the conditions for it to have N-S bifurcations. In our 
analysis, we first vary // while treating C as a constant, then vice versa, and find the value of 
H and C for which the eigenvalues have a modulus of unity and become complex conjugates. 
Then, we apply the genericity conditions for a Neimark-Sacker bifurcation to determine 
conditions on T which make the system a generic 2-D map. In section 4 we apply our 
theorem for the conditions on T to test functions that Gilet proposed. While we use the 
same form of test functions, we also study test function shapes different from those studied 
in [ ]. We end in section 5 by noting an interesting homoclinic-like bifurcation that appears 
in our simulations. 


2 Basic properties of the map 

In this section we find the fixed points as done in [;]. Then we analyze the stability of the 
fixed points for two cases: holding C constant and holding /j constant. It should be noted 
that a family of fixed points is defined as a set of multiple fixed points, and all calculations 
are done without loss of generality. 

First, define F : M 2 >->■ M 2 as, 

H[w n + '&{x n )\ 1 , , 

X n - CWn^'{x n ) _ 

We find the fixed points of (1) by solving, 


w n+ i 

% n -\-1 


F(w n} X n ) 


w * = (i[w* + T(x*)] 
x* = x* — Cw*\ I/'(a;*). 
Notice that this gives us two families of fixed points: 


(0,:r*) such that T(a;*) =0; n ^ 0 and 
(w*,x*) such that w* = —^—T(:r*). 4/'(x*) = 0; C ^ 0,/i ^ 1. 

1-p, 
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(3) 

(4) 







The derivative matrix for F is, 


DF(w,x)=( Jins i ***'(% s ) ■ (5) 

v ’ ’ \ -CS\x) 1 - Civ'S" (x) J y ’ 

2.1 Properties of F in neighborhoods of the second family 

We first notice that for fixed points (4), // e [0,1) and C € (0,1] because if /x = 1 or C — 0 
this fixed point does not exist. Substituting the fixed points (4) into the derivative matrix 
(5), we obtain 


DF(w*,x*) 

Therefore, the eigenvalues are 


/j 0 


Ai = /x, A 2 = 1 - ^ 'S(x*)'S"(x*) 

1 — /i 

Observe that since fi G [0,1), if T(£*)T"(x*) < 0 the fixed point is a saddle. However, if 
T(a;*)T // (£*) > 0, we always have A 2 < 1, so we must find a condition for which A 2 = —1. 

If we hold C constant we derive the condition /t, 

1 

^ “ l + C , T(xd^"(£*)/2' 

Then, for fi < fi this fixed point is a sink and for fi > /t this fixed point is a saddle. If // = /i, 
the fixed point becomes nonhyperbolic. 

If we hold /i constant we derive the condition (7, 


(7 = 


2(1 ~h) 
/ uT(£*)T // (a:* 


( 6 ) 


If C < C the fixed points are stable, but if C > C the fixed points are saddles. When C = C, 
A 2 = —1, which indicates the possibility of a flip bifurcation. 

When the fixed point is a saddle, it is easy to see that the linear stable and unstable 
manifolds are, 


W L = {(w,x) : x = £*}, (7) 

Wy m = {(w } x) \w = w*YZ ( 8 ) 


2.2 Properties of F in neighborhoods of the first family 

Since we can translate any fixed point (3) to the origin, without loss of generality, we assume 
the fixed point is (in*,#*) = (0,0) with the relevant conditions on Tf Substituting the fixed 
point (3) into the derivative matrix (5), we obtain 


DF{ 0,0) 


V A^'(O) \ 
-C^'(O) 1 J ' 


(9) 
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The characteristic polynomial is 

G u - A)(l - A) + C7iT'(0) 2 = 0 => A 2 - (1 + /i)A + /x(l + CT'(O) 2 ) = 0, 


and the eigenvalues are, 


A = 1(1 + V.) ± l - y4^(l + C>P'(0) 2 ) - (1 + W- (10) 

Notice that A G C \ M if 4/x(l + (TT^O) 2 ) — (1 + /x) 2 > 0. 

Holding C constant gives us the condition, 


2CT'(0) 2 - 2a/CT / (0) 2 (C'T / (0) 2 + 1) + 1 < /x < 2CT'(0) 2 + 2v'C'^ , (0) 2 (C$ , (0) 2 + 1) + 1 

( 11 ) 

Now, if | A| < 1, we get a stable focus, and when |A| > 1 we get an unstable focus. This 
indicates that there may be a N-S bifurcation at the fixed point when |A| passes through 
unity, which occurs when /x goes from /x < /x* := 1/(1 + CT^O) 2 ) to /x > /x*. 

If we hold /x constant we have, 


C > 


1 

T'(0) 2 


~ (l + h) 2 
4/i 


( 12 ) 


and plugging in this condition shows that the fixed points are always stable when the eigen¬ 
value is real. When the eigenvalues are purely complex conjugates the fixed points undergo 
N-S bifurcations, similar to the constant C case. The eigenvalue passes through unity when 

a = (i/m - i)mo) 2 . 


3 Neimark—Sacker bifurcation 

In a continuous two-dimensional dynamical system, a Hopf bifurcation occurs when a spiral 
fixed point changes stability and produces a limit cycle as a parameter is varied. The 
bifurcation is supercritical if a stable limit cycle arises as the parameter is varied forward. 
Conversely, the bifurcation is subcritical if an unstable limit cycle is created as we vary the 
parameter backward. Similarly, for a discrete two-dimensional map a N-S bifurcation occurs 
in the same manner. Due to the discrete nature of the dynamical system, an invariant closed 
Jordan curve (a topological circle) is born instead of a limit cycle. 

The existence of a N-S bifurcation at a specific fixed point for a single parameter can be 
shown by verifying the relevant genericity conditions and showing that the eigenvalues are 
complex conjugates of modulus one as the parameter passes through the critical value. The 
genericity condition is pivotal in proving a bifurcation (a topological change in the set of 
iterates) exists. Furthermore, the complex conjugate eigenvalues prove that the fixed point 
is a foci before and after the bifurcation, thereby allowing the creation of an invariant circle. 

We consider the genericity conditions outlined in [26]. That is, the map must be locally 
conjugate near the fixed point to a specific normal form, there are no strong resonances, and 
the first Lyapunov coefficient must be nonzero. The Lyapunov coefficient also determines if 
the bifurcation is supercritical (if negative) or subcritical (if positive). 

While the N-S bifurcation is a local phenomenon about a single fixed point, it should be 
noted that all calculations are done without loss of generality. Therefore, we can apply the 
theorem to individual fixed points, some of which will yield N-S bifurcations. 


6 









3.1 Neimark—Sacker bifurcation in the parameter /jl 

In [I], Gilet conjectured that a supercritical N-S bifurcation occurs at the fixed points (3). 
He also observed evidence of this in the iterates of the map for chosen test functions. Here 
we prove the map ( 2 ) is generic and a N-S bifurcation occurs as we vary the parameter /x at 
the fixed point (c.f. [26, 27, 28]). We also show the map allows for both supercritical and 
subcritical N-S bifurcations. 


Theorem 1. The map (2) is generic about some fixed point (w*,xf) if the eigenmode satisfies 
the following property, 

d =^ , "(a;*)^ , (^)(l + Ct>\x*Y){ 1 + 2CT'(o;*) 2 )(4 + 3CT'(a;*) 2 ) 

+ 2^'\x*) 2 {5 + GV Or *) 2 [1 -CT'(a ;*) 2 (31 + 21CT'(a:*) 2 )] } ^ 0. (13) 


and a N eimark-S acker bifurcation occurs at the fixed points (3) when 

1 

1 + Ctf'(x,) 2 ‘ 


(14) 


Furthermore, if d < 0, the map undergoes a supercritical Neimark-Sacker bifurcation, and if 
d > 0 , the map undergoes a subcritical Neimark-Sacker bifurcation in a neighborhood of the 
fixed point. 

Proof. As in Section 2 . 2 , without loss of generality we translate the fixed point to (w*, x*) = 
(0,0). We showed in Section 2.2 that the pair of eigenvalues A are complex conjugates if 


2C'T'(0) 2 - 2ycr(0) 2 (C'I'(0) 2 + 1) + 1 < n < 2C>t'(0) 2 + 27Cr(0) 2 (C'I"(0 ) 2 + 1) + 1 

(15) 


and | A | = 1 when 




1 

1 + CT'(0) 2 ' 


This shows that a Neimark-Sacker bifurcation occurs at the fixed point if the map is generic. 

Next we show the map is generic (c.f. [26, 27, 28]) via three conditions (C.l), (C.2), and 
(C.3), 


(C.l) We show that r'(/x*) 7 ^ 0, where r = |A| = \Jp(l + C'T'(O) 2 ) 
Notice, since r(/x*) = 1, 


i ( r(#,)2 ) 


= 2 r(/x*)r'(/x*) = 2 r'(/x*), 




so if ± (r(/x) 2 ) 


7 ^ 0, 7 ^ 0. Then, since C > 0 and ^'(0) e M, 






= (i + ct'(o) 2 ) + 0 


( 16 ) 


M=/U* 


This shows that the transversality condition is satisfied. 
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(C.2) We show the arguments of the eigenvalues satisfy the first set of nondegeneracy con¬ 
ditions for a N-S bifurcation. 

Let d* = tan -1 A, where 


^ _ Vj4/i*(l + CH/^O) 2 ) — (1 + /i*) 2 _ a/4 — (1 + fj,*) 2 
1 + /i* 1 + /i* 

Observe that 0* = 0 if A = 0, d* = ±7r if A = 0, 0* = ±27r/3 if A = ±/3, and 
d* — > ±7t/2 as A —> ±oo. Since A is clearly positive and bounded, this rules out each 
case except A = \/3. In order to get /3 we need 1 + /x — 1; however, since /x is positive 
this is not possible. 

Thus the first nondegeneracy condition is satisfied. 

(C.3) We compute the normal form for the N-S bifurcation and derive the conditions for 
which the second nondegeneracy condition is satisfied. 

First, we compute the eigenvectors: ( DF)q = A q and ( DF) T p = A p, 


Q = 


A — /x* 


and p 


For the normalization, we take 


A-! \ 

M.'P'(O) ) 


< p,q >=p-q = !/<.*'(0)V3 - (1 + /(.)- 


(18) 


(19) 


In order to simplify the calculations for the normal form let us change variables from 
M 2 to C 


w 

X 


= zq + zq = 
Substituting this into F yields 


(z + z)p^'{ 0 ) 

(A — pA)z + (A — pA>z 


( 20 ) 


F = 


; where 


F\ 

F 2 __ 

F\ ={z + ^)/x;T / (0) + /x*T((A — p*)z + (A — /x*)z), 

F -2 =(A — pA)z + (A — p*)z — C(z + 5)/x*'L , (0)T'((A — pA)z + (A — n*)z). 


( 21 ) 


Now we take the inner product (in order to transform the map into the normal form) 


< p,F >=(z + z)pl( A - l)T'(O) + ,u*(A - 1)T((A - p*)z + (A - /x*)z) + (A - /x*)/x*'L , (0 )z 
+ (A — (0)z — C (0) 2 (z + ^)T , ((A — p*)z + (A — /x*)z). (22) 



Finally, to get the normal form we divide this through by < p,q > for the sake of 
normalization and take the Taylor series in order to ignore higher order terms, 


H = 


(z + z)nl(X — l)T / (0) + /x*(A — 1) Y ~^ d z iz k ^{0)z J z k + (A - /i*)/r*'F , (0)z 


j+fc>i 


+ (A - )z -C^’(ti) 2 (z + z) ( T'(0) + Y —d zj - zk ^{0)z j z 


j-\-k> 1 


0)^4 - (1 + /'. ]- ) , 


(23) 


where 9 23 yA'T(0) := d z j z k^>((\ — n*)z + (A — /i*)£)|(z,z)=(o,o)> and similarly for d z j z k^'( 0). 
By matching linear terms, it is easy to show the normal form can be written as, 


H = Xz+ 


(A-1) 5T -c^'(a) 2 (z + z) V -Ld^wfp)z’z* 


j\k\ 

j+k> 2 J 

i*'( 0)^4-(l + /i*) 2 


j-\-k> 1 


(24) 


We are ready now to compute the nondegeneracy condition required to satisfy the final 
genericity condition. From [26, 28], we have the formula, 

d( 0) = Re + Re ®>9n) - |l3li| 2 - jM 2 . (25) 

where gjk/(j^kl) is the coefficient of the z J z k term. We compute the relevant terms 


(A - 1)(A - /i*) 2 T"(0) - 2C7r*T'(0) 2 (A - ^)T'(O) 

920 ~ i^O^T 2 ^(o) 

(A - 1)(A - /i*) 2 T"(0) - 2C/i*T / (0) 2 (A - ^)T'(O) 

902 ~ i^O^T 2 ^(o) 

_ (A — 1)(A — /i*)(A — / u*)'h"(0) C/i*'h / (0) 2 (A — /i* + A — /i*)'F / (0) 

911 ~ ?v / 4-(l + /i*) 2 ^(0) V 4 -(l + ^) 2 ^(0) 

(A - 1)(A - /i*) 2 (A - /r*)tf"'(0) Cn^'{ 0) 2 ((A - ^) 2 + 2(A - /r*)(A - /i*))T"'(0) 

(791 =- — - — - — - 

i\J 4 - (1 + /i*) 2 T / (0) i^/4 - (1 + ^( 0 ) 

If we factor out certain terms this formula greatly simplifies to 
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^"( o ) 2 Re ( iry ^ n ) ^"( 0)2 I |^ n |2 + 1|^ 02 |2 
T'(0) 2 (4 — (1 + yU*) 2 ) T'(0) 2 (4- (1 + yU*) 2 ) ’ 

(26) 


§20 — (A — //*) [A 2 + //.* — 2C4' / (0) 2 /i* — A(1 + /x*)] (27) 

g 0 2 — (A - //*) [l -//. - 2Ctf , (0) V - A/t* - A] (28) 

<7n = —1 + A+A/x* — C'4' , (0) 2 A / u* — A^/i* + A/i* — C f 4' , (0) 2 A/i* 

*- ft* — l§t + 2C , 4/ / (0) 2 /x 2 + A/x 2 (29) 

§21 = (/^* — A) [l — A — A/x* + C4/ / (0) 2 A/x* + A^/x* — A/x* + 2C f \E' , (0)^A/x* 

+/i* + /x* — SC'I'^O) 2 /! 2 — A/x 2 ] . (30) 


Then, substituting in for A and /x* gives 


d(0) = C 2 T'(0) 2 T w (0)T'(0)(1 + CT'(0) 2 )(1 + 2CT'(0) 2 )(4 + 3CT'(0) 2 ) 

+ 2T"(0) 2 {5 + CT'(O) 2 [lCT'(O) 2 (31 + 21CT'(0) 2 )] } 

/ [4(1 + CT'(0) 2 ) 4 (4 + 3CT'(0) 2 )] . (31) 

Since 

C 2 ^\ 0) 2 

4(1 + CT'(0) 2 ) 4 (4 + 3CT'(0) 2 ) > ’ 
we need only be concerned with 

d :=V'"{x*)V'{x*)(l + CT'(o;*) 2 )(l + 2CT'(a;*) 2 )(4 + 3C^'(x*) 2 ) 

+ 2T"(o;*) 2 (5 + CT'(x;*) 2 [l - CT'(o;*) 2 (31 + 21CT'(x;*) 2 )] } (32) 

Consequently, the N-S criteria imply that the bifurcation occurs if d ^ 0, and it is 
supercritical for d < 0 and subcritical for d > 0. 

This shows that for certain properties of T the second nondegeneracy condition is 
satisfied, thereby completing the proof. 

□ 
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3.2 Neimark—Sacker bifurcation in the parameter C 


In this section we prove the map (2) is generic and a N-S bifurcation occurs as we vary the 
parameter C at the fixed point. 


Corollary 1. The map (2) is generic about some fixed, point (in*, x*) if the eigenmode satisfies 
the following property, 

d = tf'(x*)tf'"(x*) (6 — fi — /i 2 ) + 2T"(x*) 2 (32/i - 6/i 3 - 21) 0. (33) 

and a N-S bifurcation occurs at the fixed points (3) when 


C = C* 



(34) 


Furthermore, if d < 0, the map undergoes a supercritical N-S bifurcation, and if d > 0, the 
map undergoes a subcritical N-S bifurcation in a neighborhood of the fixed point. 


Proof. As in Section 2.2, without loss of generality we translate the fixed point to (in*, x*) = 
(0,0). We showed in Section 2.2 that the pair of eigenvalues A are complex conjugates if 


C > 


T'(0) s 


(i + nY 


4/i 


-1 


(35) 


and | A | = 1 when 


c = a = 


i 

- - i 

h 


T'(0) 2 

This shows that a Neimark-Sacker bifurcation occurs at the fixed point if the map is generic. 

Most of the arguments for genericity follow directly from the proof of Theorem 1. For 
the first condition we use the same r and study r(C). Just as with r(/i), it suffices to study 
r(C) 2 . Then, 




= 1 (1 if <I>'(0) # 0. 


(36) 


c=c« 


Hence, if T'(0) Y 0; the map satisfies the transversality condition. 

The argument to show no strong resonances is the same as that of Theorem 1. 

For the first Lyapunov coefficient, the calculations are the same up to (30), by replacing 
C with C* and /i* with /i. Then, substituting in for A and C* gives 

(d-l? 


d( 0) = 


4T'(0) 2 (3 + /i) 


[^'(x^T'^x*) (6 - /i - /i 2 ) + 2T"(x*) 2 (32/i - 6/i 3 - 2l)] . (37) 


Since 


(/* - l ) 2 


> 0 


4T / (0) 2 (3 + /i) 

we need only be concerned with, 

d := T / (x*)T /// (x !t ,) (6 — /i — /i 2 ) + 2T"(x*) 2 (32/i — 6/i 3 — 2l) 


(38) 


Consequently, the Neimark-Sacker criteria imply that the bifurcation occurs if d Y 0, and it 
is supercritical for d < 0 and subcritical for d > 0. 

This shows that for certain properties of T the genericity conditions are satisfied, thereby 
completing the proof. 

IS 
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4 Application of theoretical results to test functions 

In [1], Gilet uses the test functions, 

T (x, (3) = —= cos (3 sin 3x H—— sin (3 sin 5x (39) 

a/TT A/7T 

where /? E [0,7r] is a fixed parameter that can be changed in order to tweak the shape of the 
eigenmode. He shows numerical constructions of the iterate-space for f3 = tt/3 and f3 = 7t/6, 
and studies the statistics of the iterates for other values of f3, all the while holding C = 0.05. 
We first follow in his footsteps and hold C = 0.05 and vary / u, then we hold /x = 0.5 and 
vary C. 

To illustrate how small the parameter regime for a subcritical N-S bifurcation is, let us 
consider d for the origin, since it is always a fixed point. We plot over (3 E [0,7r], first for /x* 
(i.e. holding C = 0.05), then for C* (i.e. holding /i = 0.5). This is shown in Fig. 2. Recall 
that d < 0 denotes a supercritical N-S and d > 0 denotes subcritical. 

d d 




Figure 2: Plot of first Lyapunov coefficient for Neimark-Sacker bifurcations in /j and C 
respectively. 


4.1 Simulation results for Neimark-Sacker bifurcation in fi 

As in [I], we first study the eigenmode for f3 = tt/3, for which the map exhibits supercritical 
N-S bifurcations at various fixed points, but we also study (3 = 57t/ 6, for which the map 
exhibits a subcritical N-S bifurcation at the origin. We hold C = 0.05, which corresponds 
to a droplet much smaller than the cavity being created by an impact, and we vary /x, which 
corresponds to the system approaching high memory. 

It should be noted that for the sake of plotting, we consider the map (2) to be written as 


71+1 

Wn+1 


F(x n , w n ) 


x n - Cw n ty'(x n ) 

fi[w n + 'F+ n )] 


(40) 
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4.1.1 Supercritical Neimark—Sacker bifurcations for (3 = 7t/3 

For /3 = 7t/ 3, we now consider the map (40) on the domain x G [0, 7t/ 2]. On this domain, the 
map has three fixed points (x*,0) about which N-S bifurcations may occur; one of which is 
the origin with the other two being the following: 


x* 




( 

15 + a /3 — \J 6(8 — \/3) ^ 

[\ 

9 — a /3 + \J 6(8 — \/3) } 

/ 

15 + y/3 + \J 6(8 — \/3) ' 

[\ 

9 — a /3 — \J 6(8 — v^3) y 


« 0.7269 


« 1.3515 


At the fixed points a supercritical N-S bifurcation occurs for //* = 0.64894, 0.742027, 
0.879451, respectively, for which d = —4076.61, —1747.52, —410.779. We also have three 
other fixed points: (^ 0.3490,^ ^^(0.3490)), (» 1.0128, « ^^(1.0128)), (tt/2, ^^(tt/2)), 
which are saddles (the w position varying with /i) and therefore not subject to N-S bifurca¬ 
tions. Let us vary the parameter /j from /j = 0.5 to /r = 0.89. We illustrate the progression 
of the bifurcations in Fig. 3. In Fig. 3a, the relevant fixed points are all stable foci. Next, 
in Fig. 3b, we have passed the critical value for the origin and a stable invariant circle is 
now visible. Finally, in Fig. 3c and 3d, the critical values for the next two fixed points, 
respectively, are passed. We notice in Fig. 3c and 3d that the preceding fixed point (s), 
respectively, are now unstable focus(foci). This shows that each fixed point undergoes a 
supercritical N-S bifurcation in accordance with Theorem 1. It should be noted that since 
the test functions are periodic, on the entire domain of x 6 1, N-S bifurcations occur at 
(countably) infinitely many fixed points with the final bifurcation occurring for fixed points 
corresponding to x* ~ 1.3515. 

Physically, the droplet about any fixed point is initially (Fig. 3a) oscillating with de¬ 
caying radius. Then, due to the nature of the single eigenmode used by Gilet, the droplet 
corresponding to a position near the origin will undergo a N-S bifurcation (Fig. 3b). This 
corresponds to the droplet oscillating about the origin with a fixed radius. Then as /j is 
varied further, the invariant circle about the origin disintegrates and it becomes unstable, 
which corresponds to the orbit of the droplet increasing radially until it reaches the neigh¬ 
borhood of one of the neighboring fixed points. Now, since a N-S bifurcation has already 
occurred at the next fixed point (Fig. 3c), the iterates will bypass this region and converge 
to the following fixed point, which corresponds to the droplet jumping to an orbit around a 
different center. Similar behavior occurs at every other N-S fixed point. 

In addition to N-S bifurcations, we notice a curious phenomenon. As the saddle fixed 
points move away from the N-S fixed points the trajectories from the neighborhood of one 
N-S fixed point crosses into that of another. Furthermore, as the invariant circle increases 
in radius it seems to collide (more precisely, gets arbitrarily close) with a stable manifold of 
one of its neighboring saddles just before the onset of chaos. We shall analyze this further 
in Section 5. 
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(c) Invariant circle about (~ . 726 ! : fj, = .79 


(d) Invariant circle about (« 1.3515,0): /u = .89 


Figure 3: For each plot the green lines represent the linear stable manifold and the red 
lines represent the linear unstable manifold at the respective saddle fixed points. The black 
markers represent iterates originating from a neighborhood of the origin, the cyan markers 
represent iterates originating from a neighborhood of ya .7269 0), and the magenta markers 
represent the iterates originating from a neighborhood of (~ 1.3515,0). 
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4.1.2 Subcritical Neimark—Sacker bifurcations for (3 = 57t/6 

For /3 = 57 t/ 6 let us once again consider the map (40), however let us simplify matters and 
restrict our domain to a neighborhood of the origin. Here lies a subcritical N-S bifurcation 
at /x* = 0.999847 where d = 4.88756. Since this is subcritical, let us vary our parameter // 
backwards from /x = 0.9999 to fi = .999. This is illustrated in Fig. 4. First the origin is an 
unstable focus, then as /x passes /x* backwards the origin becomes a stable focus, however if 
an initial point is taken further out the iterates diverge, which indicates an unstable invariant 
circle. In Fig. 4b, we represent this unstable invariant circle by the black dashed curve. 



(a) Unstable focus: /x = 0.9999 (b) Unstable invariant circle: /x = 0.999 

Figure 4: A star denotes the initial point. In the second plot the blue iterates have an 
initial point inside the invariant circle, and the red iterates have an initial point outside the 
invariant circle. Finally, the black dashed curve represents the unstable invariant circle. 

While from a dynamical systems perspective, we vary the parameter backward for a 
subcritical bifurcation, this would not make sense for physical intuition. Rather let us begin 
at Fig. 4b, where a droplet close enough to the origin would oscillate with a decaying radius, 
however if the droplet were exactly on the unstable invariant circle it would remain at that 
constant radius, and if the droplet were slightly further out it would oscillate with a growing 
radius and escape the neighborhood of the origin. Now as we increase the memory the 
invariant circle shrinks until the critical parameter value /x* = 0.999847, when it ceases to 
exist and the origin becomes unstable, which corresponds to the droplet anywhere near the 
origin to oscillate with increasing radius. It should be noted that the growth of the radius is 
not indefinite, and at some point it may be drawn to another fixed point or invariant circle, 
which is usually the case. 
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4.2 Simulation results for Neimark—Sacker bifurcation in C 

Varying C produces quite exotic trajectories, one example of which we present in this section, 
holding 11 = 0.5 and f3 = 7t/ 3. Consider two fixed points, (x*, 0) such that 


/ 


x* = — tan 


-l 


15 + v ^3 + yj 6(8 — \/ 3 ) ^ 
^ 9 — a /3 — \ J 6(8 — a / 3 ) } 


-1.3515, 


/ 


x* = —7T + tan 


-l 


15 + v ^- \ j 6(8 — \/ 3 ) ^ 
9 — a /3 + \ J 6(8 — \/ 3 ) y 


-1.7901. 


At these fixed points a supercritical Neimark-Sacker bifurcations occurs simultaneously 
due to the symmetry for 77* = 0.36477, where d = —578.487. The other fixed point is 
(—7 t/2, T(-7t/ 2)), which is a saddle. 

For this case, varying C forwards corresponds to the drop size increasing, which leads 
to greater inertia. With increased inertia, the distance the droplet travels due to an impact 
is greater, and hence sensitive dependence may come into play. First, the N-S bifurcation 
is illustrated in Fig. 5. Here the physical interpretation of the bifurcation is the same as 
in section 4.1, the only difference being the change in the size of the droplet as opposed to 
changes in memory. 



(a) Stable foci: C = 0.35 (b) Stable invariant circle: C = 0.4 

Figure 5: The red trajectories start on the left of the saddle blue trajectories start on the 
right of the saddle. 


Now, as we increase C, we obtain the exotic trajectories as illustrated in Fig. 6. We 
observe the characteristic stretching and folding of a chaotic system, which may lead to 
fractal structures. The way in which the invariant circles approach the stable manifold that 
separates them shows evidence of what appears to be an unusual and possibly new type of 
global bifurcation, which is briefly discussed further in the next section. 
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(a) Deformed invariant circles: C = 0.5 (b) Evidence of global bifurcation: C = 0.6 

Figure 6: The red trajectories start on the left of the saddle blue trajectories start on the 
right of the saddle. 

5 Evidence of global bifurcations leading to chaos 

In Fig. 3d we observe the appearance of quite exotic dynamics. If we continue to vary p 
we find more exotic, chaotic-like, dynamics. We illustrate this in Fig. 7 where we see a 
scatter of iterates in a manner indicative of chaos. Moreover, we observe as p is varied, the 
iterates remain within a compact set. Here the iterates seem to satisfy the transitivity and 
sensitivity conditions for chaos. So far, all we have are simulations and informed intuition: 
further investigation is clearly required in order to obtain a more precise characterization 
these dynamical properties. 

Now let us vary p near the onset of chaos to study the bifurcation. This is illustrated in 
Fig. 8. First (Fig. 8a) we observe our invariant circles for p = 0.913 as before, then as we 
increase to p = 0.914 the invariant circles start to disintegrate (Fig. 8b). Finally, in Fig. 8c, 
for p = 0.915 we see the iterates bouncing between the neighborhoods of their respective N-S 
fixed point and perhaps nested invariant circles (both stable and unstable). This provides 
evidence of an exotic new bifurcation as the invariant circle approaches the stable manifold 
of the saddle and perhaps even a cascading N-S doubling bifurcations, which is an infinite 
succession of stability changes of invariant circles that create new pairs of invariant circles 
(see, e.g. [29], where they are called cascading Hopf doubling bifurcations). 

Furthermore, we notice in our figures that iterates prefer to intersect the linear stable 
manifolds. If we do not restrict the region of the plot, we observe that the iterates seem 
to prefer certain positions as shown in Fig. 9. In [l], Gilet gives a detailed study of the 
statistics associated with the positions of the iterates and in the simulations he observes a 
probability density function inversely proportional to |T / (x)|. While we do not study the 
statistics, we observe that the linear stable manifolds pass through the roots of T^a:), which 
hints at connections between the dynamical systems aspects and statistical analysis of this 
model. 
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(a) Evidence of chaos: fi = 0.94 (b) Evidence of chaos: n = 0.96 


Figure 7: Exotic, chaotic-like, trajectories. For each plot the i,reer lines represent the linear 
stable manifold and the red lines represent the linear unstable manifold at the respective 
saddle fixed points. The black markers represent iterates originating from a neighborhood 
of the origin, the cy; markers represent iterates originating from a neighborhood of 

0), and the magenta markers represent the iterates originating from a neighborhood 
of (« 1.3515,0). 
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(a) Invariant circle: fj, = 0.913 (b) Invariant circle disintegrat- (c) Iterates after collision: p = 

ing: /r = 0.914 0.915 

Figure 8: The green lines represent the linear stable manifoh and the red lines represent 
the linear unstable manifold at the respective saddle fixed points. The blue markers are the 
iterates. 
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Figure 9: Iterates of the map (40) generated from four initial points chosen near the origin. 


6 Conclusions 

Previous investigations of the walking phenomena have shown evidence of exotic dynamics, 
such as Hopf bifurcations and chaos, mainly through experiments and numerical simula¬ 
tions. However, the equations have generally been too complex to prove the existence of 
certain bifurcations, chaos, and other topological properties. In recent years, discrete mod¬ 
els have been studied by Fort et al. [12], Eddi et al. [17], Shirokoff [23], and Gilet [1]. The 
models of Shirokoff [23] and Gilet [ ] are discrete dynamical systems, which are often more 
accessible to rigorous analysis. In this paper we prove some of the properties observed and 
conjectured for Gilet’s model; namely, the existence of N-S bifurcations. Similar dynamical 
systems analysis is likely in future to prove additional results concerning bifurcations, chaos, 
structural stability, and other properties. Proving these properties is not only important 
mathematically, but may also provide insight into the phenomena observed for walkers and 
for their continuous models. 

We proved Gilct’s conjecture [ ] about the existence of supercritical N-S bifurcations 
when varying the parameter /j. In doing so, we are also able to prove the existence of 
subcritical N-S bifurcations in the parameter /i, which had been missed by the original 
simulations. We then varied the parameter C, which had not been studied before, and 
proved the existence of N-S bifurcations there as well. 

In order to verify the validity of our theory, the model was simulated using test functions 
for T. For the sake of consistency the test functions used were of the same form as [I], but the 
shapes were changed by varying (3. By doing so were able to analyze the transition between 
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supercritical and subcritical N-S bifurcations (Fig. 2) for both /j and C. We reproduced 
one of the simulations done in [1] in order to show consistency between the two works. Then 
we ran additional simulations of the various cases in our theorems, all of which demonstrate 
complete agreement between our results and the numerics. 

In addition, we observed evidence of more exotic phenomena, which seem to lead to 
chaos. One such curious phenomenon is a homoclinic-like bifurcation when the invariant 
circles get ever closer to the stable manifold of a neighboring saddle fixed point. This also 
appears to give rise to other invariant circles leading to cascading N-S bifurcations wherein 
the invariant circles change stability and give birth to twin stable invariant circles. 

Finally, while not the focus of this work, Gilet, as did we, observed interesting statistics for 
the position of the iterates. We notice that after the onset of chaos the iterates prefer certain 
positions, namely along the stable manifolds of each saddle fixed point. In [1], Gilet studies 
the statistics in detail through the use of approximate probability density functions. It may 
be possible to obtain a more analytical characterization using measurable dynamics theory, 
which is a problem worthy of future investigation. He writes, “Future work could include 
the search for wave-particle coupling dynamics that yields PDF (a;) ~ ^(a;)! 2 as observed in 
quantum mechanics.” We agree that this would be a fruitful endeavor and perhaps studying 
the dynamics of the system further will provide insight into the right choice of wave-particle 
coupling. 
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